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We report the observation of e^e~ <f>r] near y/s = 10.58 GeV with 6.5 a significance in the 
K'^ K~^^ final state in a data sample of 224 fb~^ collected by the BABAR experiment at the PEP- 
II e^e~ storage rings. We measure the restricted radiation-corrected cross section to be cj{e^e~^ 
(j>ri) =2.1 ± 0.4(stat) ± O.l(syst) fb within the range | cos^*| < 0.8, where 9* is the center-of-mass 
polar angle of the 4> meson. The (j) meson is required to be in the invariant mass range of 1.008 
< < 1.035 GeV/c^. The radiation-corrected cross section in the full cos 9* range is extrapolated 
to be 2.9 ± 0.5(stat) ± O.l(syst) fb. 

PACS numbers: 13.66.Bc, 13.25.-k, 14.40.Ev 



The large data samples of the B factories provide 
an opportunity to explore rare exclusive quasi-two-body 
processes in e^e~ annihilation, such as final states pro- 
duced through one virtual photon with negative C-parity 
{J/iprjc or other double charmonium states) [l], and 
two- virtual-photon annihilation (TVPA) with positive C- 
parity (p°p° or The process e+e^ J 0^110 and 

other double charmonium processes are observed at rates 
approximately ten times larger than the expectation from 
QCD-based models Various theoretical efforts have 
been made to resolve the discrepancy between experi- 
mental and theoretical results [1, [EI3- Another avenue 
to explore this puzzle is provided by the related process 
e+e" (j)rj. A recent observation of '0(3770) (f>r] at 
a branching fraction of (3.1 ± 0.6 ± 0.3 ± 0.1) x 10"^ ^ 
also stimulates a search for T{4S) (prj. 

We report the observation of e^e^ — )■ (jji], which is 
analogous, in the s quark sector, to the process e'^e" 
J/V'??C7 since the rj meson has an ss quark-pair com- 
ponent. The Feynman diagram for the most likely pro- 
duction mechanism is shown in Fig. [TJ However, since 
rj is not purely ss, the cross section for this production 
mechanism is determined by the projection onto the ss 
component of the rj meson. A calculation using the QCD- 
based light cone method with relativistic treatment for 
the light s quark is possible and therefore can provide a 
theoretical estimation [Oj. 

The (j>ri combination is a vector — pseudoscalar (VP) 
final state. The production rates for e+e^ VP can 
be described by form factors, which are predicted in 
QCD-based models [13, [HI, [I2|- Different models pre- 
dict different dependences on center-of-mass (CM) en- 
ergy squared s. The recent measurements of e~^e 



VP [ujTT^ T pr] and prj') from BES 13|, investigated the 
s dependence of the cross sections and form factors in the 
energy range from 3.65 to 3.773 GeV. It is interesting to 
investigate the s dependence over a wider energy range. 
Since CLEO measured the cross section for e+e^ — > ^77 
at CM energy ^/s = 3.67 GeV [1], a measurement of the 
same process at ^/s — 10.58 GeV provides a meaningful 
test of the s dependence. 

This analysis uses 204 fb^^ of e+e^ colliding beam 
data collected on the T(4S') resonance at ^/s = 10.58 
GeV and 20 fb"! collected 40MeV below the r(4S) 
mass with the BABm detector at the SLAG PEP-II 
asymmetric-energy B factory. The BaBAR detector is 



FIG. 1: Possible Feynman diagram for e^e (j)rj. 



described in detail elsewhere [15|. Charged-particle mo- 
menta and energy loss are measured in the tracking sys- 
tem that consists of a silicon vertex tracker (SVT) and 
a drift chamber (DCH). Electrons and photons are de- 
tected in a CsI(Tl) calorimeter (EMC). An internally 
reflecting ring-imaging Cherenkov detector (DIRC) pro- 
vides charged particle identification (PID). An instru- 
mented magnetic flux return (IFR) provides identifica- 
tion of muons. Kaon and pion candidates are identi- 
fied using likelihoods of particle hypotheses calculated 
from the specific ionization in the DCH and SVT and 
the Cherenkov angle measured in the DIRC. Photons are 
identified by shower shape and lack of associated tracks. 

To reconstruct ^77 in the K~^j mode, events with 
exactly two well-reconstructed, oppositely charged tracks 
and at least two well-identified photons are selected. 
Charged tracks are required to have at least 12 DCH 
hits and a laboratory polar angle within the SVT accep- 
tance, 0.41 < 6 < 2.54 radians. The laboratory momenta 
of the kaon candidates are required to be greater than 
800 MeV/c to reduce background. The two tracks se- 
lected must both be identified as kaons. We fit the two 
tracks to a common vertex, and require the probabil- 
ity to exceed 0.1%. The photon candidates are required 
to have a minimum laboratory energy of 500 MeV. The 
invariant mass distribution of K^"fj^ after requir- 
ing the invariant mass of KK to be near the </> mass 
{rriKK < 1-1 GeV/c^) and that of 77 to be near the 77 
mass (0.4 < m^^ < 0.8 GeV/c^) is shown in Fig. [H (a). 
We accept events with a reconstructed invariant mass of 
K+R-^-f within 230 MeV/c^ of the e+e" CM energy. 
There is at most one entry per event in the region of 
interest. 

Figure [2Ib) shows the scatter plot of invariant masses 
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FIG. 2: (a) Distribution of the invariant mass (T(4S') data) 
for the K~'y^ final state near the cfirj region. The accepted 
signal region is indicated by the lines, (b) Scatter plot of the 
invariant masses of the K~ and 77 pairs for those events 
in the accepted signal region. 
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FIG. 3: Mass projections for (a) K^K pairs and (b) 
pairs in K~jj events. 
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of K^K~ and 77 pairs from the accepted e+e" 
K^K~jj events. The concentration of events indicates 
(j)ri production. 

We use a two-dimensional log-likelihood fit to extract 
the signal for the reaction e~^e~ (f/rj. Due to the fact 
that the final state particle masses are far below the e+e~ 
collision energy, we may treat the two-body masses as 
uncorrelated. Justified by Fig. [21(b), the signal probabil- 
ity density function (PDF) is constructed as a product 
of two one-dimensional PDFs, one for each resonance. 
We use a P-wave rclativistic Breit-Wigner formula to 
construct a PDF for the (f> resonance and a Gaussian 
function to model the 77 resonance. A threshold function 
q'^ /{I + q^Rt) is used to model the background in the 
system, where q is the daughter momentum in 
the (j) rest frame and Rt is a shape parameter. A linear 
function (po +Pi ■ rn^^) is used to model the background 
under the 77. 

In the fit to data, we fix the mass and width of the (p 
and the mass of the rj to the world average values 
The width of the 77, 13.6 MeV, is fixed to the resolu- 
tion obtained from simulation. The floating parameters 
in the fit are: Rt, po, pi, and the numbers of events 
for all components-^?7, (f)jj and rjK^K^. The mass 
projections in KK and 77 from the two-dimensional fit 
are shown in Fig. [3] (a) and (b), respectively. We de- 
fine the (j) mass window as 1.008 < tukk < 1.035 
GeV/c^ to reduce the systematic uncertainty due to the 
long tail of (j) masses. The extracted number of (f)r] sig- 



nal events is 24 ± 5 in the (j> mass window, with 20±5 
in the on-resonance sample and 3±2 in the off-resonance 
sample. The number of background events within the 
(j) mass window and within three standard deviations 
of the 77 mass is 7±2. The significance is estimated by 
the log-likehhood difference between signal {ln{Ls)) and 
null {ln{Ln)) hypotheses (no (p-q signal component in the 
PDF), ^j2ln(Ls/Ln), which gives 6.5 standard devia- 
tions. 

Given the negative C-parity of the (prj final state, we 
assume 4>rj is produced through one- virtual-photon anni- 
hilation. The angular distributions of cpr] from a = 1^ 
initial state, in the helicity basis [3], can be calculated 
to be: 

dN 



d cos 9* d cos 



cx sin^ 6*0 ( 1 -|-cos^ 9* +cos 2(p^ sin^ 9*), 

(1) 

where the production angle 9* is defined as the angle be- 
tween the (f> meson direction and incident e~ beam in 
the CM frame. The (p helicity angle 9^ is defined as the 
polar angle, measured in the (j) rest frame, of the mo- 
mentum direction with respect to an axis that is aligned 
with the (p momentum direction in the laboratory frame. 
The variable ip^ is the K'^ azimuthal angle around the 
direction of the (j) measured with respect to the plane 
formed by the (j) and the incoming electron. The helic- 
ity and azimuthal angles of the pseudoscalar 77 are fiat 
and thus not included in equation [TJ Integrating over 
the other two angles, the distributions of the production 
angle, (p helicity and (fi azimuthal angle are expected to 
be 1 -t- cos^ 9*, sin^ and 2 -|- cos2(/?0, respectively. The 
observed angular distributions from e~^e~ — > ^77 data are 
consistent with the above expectation but the constraints 
on these angular distributions are limited by statistics. 

The systematic uncertainty from the two-dimensional 
fit is estimated from the difference in yield obtained by 
fioating the mean, width, and resolution parameters in 
the fit. The systematic uncertainties due to PID, track- 
ing, and photon efficiency are estimated based on mea- 
surements from control data samples. The possible back- 
ground from related modes with an extra ir^ was esti- 
mated to be small (< 1%) by using extrapolations from 
statistically limited four-particle mass sidebands and we 
ignore it. The systematic uncertainties are summarized 
in Table [J 

The radiation-corrected cross section for e+e^ — s- (prj 
is calculated from: 



Observed 



Lx B{p^ KK) X B{r] 77) x e^^C' x (1 (5) 

where iVobsorvod is the extracted number of (prj signal 
events from on- and off-resonance data, L is the inte- 
grated luminosity, B{(j) KK) is the branching frac- 
tion of ^ KK, B{rj 77) is the branching frac- 



tion of 77 



77, e 



MC 



is the signal efficiency obtained 
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TABLE I: Systematic uncertainties on the cross section of (j>V- 



Source 


Systematic uncertainty % 


Photon efficiency 


3.6 


Two-dimensional fit 


1.3 


Particle Identification 


3.0 


Tracking efficiency 


2.6 


Luminosity 


2.0 


Total 


6.0 



from Monte Carlo simulation (MC), and 6 is the radi- 
ation correction calculated according to Ref. 17|. We 
obtain (1 + ^) = 0.768. The uncertainties due to the 



theoretical model and the s dependence are negligible. 
The signal MC events are generated uniformly in phase 
space. For the determination of signal cross sections, the 
MC cos 9*, cos 00 and ip^ distributions are re- weighted 
using equation [TJ The signal efficiency in the fiducial re- 
gion of I cos 0*1 < 0.8 for 077 without radiative correction 
is estimated to be 34.3%, including corrections to MC 
simulation for PID and tracking. Taking the branch- 
ing fraction of ^K^K~ as 49.1%, and 77 ^ 77 as 
39.4% 0, the final radiation-corrected cross section for 
1.008 < TO0 < 1.035 GeV/c2 within |cosr| < 0.8 near 
yi= 10.58 GeV is: 

CTfid(e+e" (pif) = 2.1 ± 0.4(stat) ± O.l(syst) fb. 

The cross section within cos 0* e [—0.8, 0.8] can be scaled 
to cos 9* E [—1.0, 1.0] by assuming a 1 + cos^ 9* distribu- 
tion to obtain: 

cr(e+e" (l>T]) = 2.9 ± 0.5(stat) ± O.l(syst) fb. 

To study the possibility that the observed signal is 
due to T(4S') decay, we scale the off-resonance signal to 
the on-resonance luminosity, and subtract it from the 
on-resonance signal. The resulting number of events, 
— 10 ± 21, is consistent with zero. The corresponding 
branching fraction for T{AS) (jir] is (-0.9±1.8) x 10"^. 
Assuming this uncertainty can be treated as Gaussian 
and normalizing to the physical region (> 0), the 90% 
confidence level upper limit is 2.5 x 10~^. 

There is currently no direct prediction for the cross sec- 
tion of this process at this energy, but the e+e" VP 



cross section is expected to have 1/s^ [13] or 1/s^ 13, 11 ] 
dependence in QCD-based models. A comparison be- 
tween our result and that of CLEO, (cr = 2.ll^;2±0.2 pb) 
at = 3.67 GeV Q, favors a l/s^ dependence (Fig.H]). 
We quantify the degree to which 1/s'* scaling is disfavored 
by scaling our measured cross section in this fashion to 
i/i = 3.67 GeV, and comparing it to the CLEO measure- 
ment. Note, however, that if CLEO did have a downward 
statistical fluctuation, both their central value and their 
uncertainty would be low. Accordingly, the uncertainty 
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FIG. 4: Cross section extrapolations based on BABAR^s mea- 
surement at y/s — 10.58 GeV assuming (black) or 
(red) energy dependence. The bands show one standard de- 
viation uncertainties in the extrapolations. The CLEO mea- 
surement at y/s — 3.67 GeV is also shown. 



we use in this comparison is the CLEO one scaled by 
the square root of the ratio, 2.6, of the predicted to the 
observed cross sections. The resulting disagreement with 
1/s* scaling is approximately 2 standard deviations. 

The form of the s dependence has important theo- 
retical implications, which may affect a wide range of 
QCD-based processes such as e~^e~ VP [13], exclu- 
sive hadronic B decays [l^ , and charmonium decays [13] . 
The large initial-state radiation sample at the B factories 
can provide another route to test the s dependence over a 
wider energy range. A direct comparison of the absolute 
cross section with a possible theoretical calculation 9] is 
also interesting. 

In summary, we have observed the exclusive produc- 
tion of (j>ri in e+e^ interactions at y/s — 10.58 GeV. Com- 
bining with CLEO's measurement and interpreting our 
result as continuum production, the measured (j)ri cross 
section favors l/s"^ dependence, which is in conflict with 
some QCD-based predictions. The 90% confidence level 
upper limit on the branching fraction B{T{AS) — > (jyrj) is 
2.5 X 10-^. 
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